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The mean-field approach of glassy dynamics successfully describes systems which are out-of- 
equilibrium in their low temperature phase. In some cases an aging behaviour is found, with no 
stationary regime ever reached. In the presence of dissipative forces however, the dynamics is 
indeed stationary, but still out-of-equilibrium, as inferred by a significant violation of the fluctuation 
dissipation theorem. The mean-field dynamics of a particle in a random but short-range correlated 
environment, offers the opportunity of observing both the aging and driven stationary regimes. Using 
a geometrical approach previously introduced by the author, we study here the relation between 
these two situations, in the pure relaxational limit, i.e. the zero temperature case. In the stationary 
regime, the velocity («)- force (F) characteristics is a power law v ~ F 4 , while the characteristic 
times scale like powers of v, in agreement with an early proposal by Horner. The cross-over between 
the aging, linear-response regime and the non-linear stationary regime is smooth, and we propose a 
parametrisation of the correlation functions valid in both cases, by means of an "effective time". We 
conclude that aging and non-linear response are dual manifestations of a single out-of-equilibrium 
state, which might be a generic situation. 

PACS numbers : 05.70.Ln, 64.70.Pf, 75.10.Nr, 83.50.Gd 



I. INTRODUCTION 



Thermal equilibrium is the situation where all fast processes have already taken place while slow processes have not 
yet started happening At the opposite, systems with slow dynamics are characterised by a broad distribution of 
relaxation times, ranging from the microscopic scale (10 -12 s) to the macroscopic one (hours or days). For instance, 
glassy systems have an equilibration time, either infinite, or much longer than the laboratory time scale. These 
systems reveal their out-of-equilibrium state in phenomena like aging or non-linear response. 

In these systems, the microscopic time scale is not the only relevant one, and much slower processes also take place. 
The slow dynamics is generally attributed to the presence of thermally activated barrier crossing in the configuration 
space, but others mechanisms, such as the so-called "entropic barriers" may also contribute H. Glassy dynamics 
is observed when the relaxation time r becomes larger than the laboratory typical time scale, as it is the case for 
supercooled liquids ||. 

The aging behaviour of spin glasses has been thoroughly investigated Q. The thermoremanent magnetization of 
field cooled samples shows a strong waiting-time dependence (where the waiting-time t w is the time interval between 
the temperature quench and the measurement). These systems have an a-priori infinite internal relaxation time, and 
the late stage of the relaxation is instead controlled by the waiting-time itself. Moreover, field cooled (waiting for 
t w and switching off the field) and zero field cooled (waiting for t w and switching on the field) show a remarkable 
complementarity of the magnetization curves [^|,^). While out-of-equilibrium, as indicated by its significant waiting- 
time dependence, the response of the system is linear in the applied field, provided this one is weak enough. 

Glassy dynamics is also observed in the dissipative dynamics of high-Tc superconductors. Supraconducting samples 
with quenched disorder, at magnetic field and temperature large enough, offer a significant resistance to a flowing dc 
current, due to the thermal motion of the flux lines. A transition line is believed to separate an ohmic regime (the 
vortex liquid) from a true superconducting state (the vortex glass) In the latter, and in the limit of a vanishingly 
small current j, the dissipation occurs by activation of "bundles" of flux lines over pinning energy barriers. According 
to the scaling theory of the vortex glass, the typical time needed for such a move, r(j), diverges exponentially fast 
as j tend to zero [|[. In this situation, the response (the voltage) is a non-linear function of the driving force (the 
current). The system is out-of-equilibrium, because of a constant rate dissipation, but stationary, at variance with 
the spin glass aging. The relaxation time r(j) which would be infinite in the absence of driving force, is regularized 
by any small but constant j, and inversely related to the magnitude of j. 
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Important and related issues are the supercooled liquids dynamics and the rheology of soft glassy materials. In the 
hrst case, the mode-coupling approach predicts an increase of the structural relaxation time r upon cooling Aging 
have been found during the early stage of molecular dynamics studies of the Lennard- Jones fluid |l(]Jll| . On the other 
hand, a constant shear rate flow seems to be able to change the value of r, resulting in a shear rate dependent viscosity, 
e.g. a non-linear response of the fluid @. More generally, this shear-thinning behaviour is well known in the context 
of soft- matter rheology, and a phenomenological treatment of this phenomenon based upon glassy dynamics has been 
proposed [ p^|Jl4| . These are situations where aging and non-linear response certainly coexist as manifestations of a 
more general glassy dynamics. 

Among the existing theoretical approaches on glassy systems, the mean-field dynamics is a very promising one. It 
was already used to suggest that the presence of dissipative forces generically prevents the aging phenomenon |]l5| , |l6| . 
In this framework, the out-of-equilibrium character of a system is made precise by the existence of a generalised 
fluctuation dissipation theorem, related in turn to the entropy creation rate p7| . 

A model of particular interest is the mean-field dynamics of a particle with a quenched pinning potential. Isolated, 
the particle presents an aging behaviour with a logarithmic growth of the time correlation functions Q,^9|. In the 
presence of a time-independent driving force, the dynamics is believed to be stationary, with a power law dependence 
of the particle's velocity in the applied force |20|. 

Recently, the author presented a geometrical description of the aging and linear response regime of this model, at 
zero temperature This approach is extended, in this paper, to the non-linear stationary regime. As a result, we 
find that aging, linear response dynamics on the one hand, and stationary, non-linear response dynamics on the other 
hand, are indeed dual manifestations of a single out-of-equilibrium state. The constant force is found to interrupt 
efficiently the aging relaxational dynamics, and to control the characteristic times, which in turn control the effective 
friction coefficient in the stationary regime. The resulting velocity-force characteristics is v oc F 4 , while the cross-over 
time between aging and stationary regime is tp cx F~ 3 . These predictions are confronted to the numerical integration 
of the mean-field equations. We finally suggest a scaling behaviour for the correlation function of this model which, 
according to our numerical findings, interpolate smoothly between the two different regimes, demonstrating their 
common origin. 



II. MEAN FIELD EQUATIONS AND THE HORNER RESULT 

We focus on the zero temperature relaxational dynamics of a particle in a quenched random gaussian potential . 
The particle evolves in a TV-dimensional space, under the simultaneous effect of a pinning force -V^ and a constant 
force and the equation of motion for the vector position x(t) is : 

x(t) = -V7(x(t)) +J 7 . (2.1) 

The potential V(x) is a quenched disorder, chosen from a gaussian distribution, with correlations (the overline stands 
for the average over the quenched disorder) : 



V(x) • V(x') = N ■ exp - " J ; V(x) = 0. (2.2) 

This form ensures a meaningful N — > oo limit, in which each coordinate Xj(£), or gradient component 9jV(x), remains 
of order one, while the norms ||x(t)||, ||VV|| scale like N 1 / 2 . The force is directed along the direction 1. 

The thermodynamic limit N — > oo is taken first, which makes the zero temperature dynamics non trivial |2^ ]. In 
this mean-field limit, the relaxation process is completely described by the displacement u, the response function r 
and the correlation functions b and d, ix being the Martin-Siggia-Rose auxiliary time p3| , pi[ . 

u(t) = N- 1 ' 2 xi(t); (2.3) 

N 

r(t,t')=N- l Y / x j (t)-ix j (t'); (2.4) 

3=1 
N 

b(t, o = n- 1 J2 fo(*)-^-(t')) 2 ; (2-5) 

N 
3=1 

= b{t,t') + [u{t)-u(t')} 2 . (2.6) 
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The Dyson equations for r, b, d, u are a closed system of coupled integro-differential equations, which contains the 
equations of [Q^o) as a particular case : 

d t r(t,t') = 6(t-t') 

-4 f ds exp(-d(M)) r(t,s) [r(t,t') - r(s,t% (2.7) 
Jo 

d t b{t,t') = (2T)-4 / ds exp(-d(t,«)) [r(t, s) - r(f', a)] 
Jo 

-4 / ds exp(-d(i,s)) r(t,s) [6(i,s) + 6(t,t') - b(s,t')]; (2.8) 
Jo 

9tu(t) = F - 4 / ds exp(-d(t, s)) r(i, s) • - w(s)]. 

(2.9) 
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The temperature term (2T) is actually zero in our case. It is also convenient to define the integrated response 1Z and 
the energy £: 

n{t,t')= J^dsr(t,s); (2.10) 

£(t) = -2 [ ds exp(-d(t,s))r(t,s), (2.11) 
Jo 

In a seminal paper, Horner described the stationary state reached by the system when driven with a finite force F, 



in the case of short range, power-law correlations ]20J]. We have found that the system (2.9) does indeed lead to 
a stationary situation, that we study in detail in this paper. One must mention however that the stationary state 
reached by the particle depends on how the system is prepared, in the same way as thermalised initial conditions can 
prevent aging in the p-spin case p5| . The stationary state is found only if the system is quen ched from a high enough 
temperature [^6|. Results for a similar driven system have also been recently published 27 



Let us summarise the main properties of the stationary solution found by Horner pC|, The correlation functions 
are time-translationally invariant (TT I) : r(t, t') — R{t — t'), B(t, t') = B(t — t') while the displacement goes linearly 
with time u(t) — v ■ t. The system ( |2.9D becomes a set of non-causal equations to be solved self-consistently. A 
non trivial feature of this solution is the emergence of characteristic time scales dependent on the velocity. With the 
notations of |2C| ], t p (v) is the characteristic time for breaking the fluctuation dissipation theorem, while t a (v) controls 
the main "a" relaxation of the correlation function B{t). These characteristic times play a very similar role than the 
time scales tf and tb respectively, introduced in pjj], and we identify subsequently tf = t p , tb = t a . In the long time 
regime, t ~ i&: 



m=q + B {W))> (2 - 12) 

with, 

t f ( v ) ~ vb-W, (2.13) 

t b (v)^v^-\ 0<r?,C< 1. (2.14) 



The exponents depend (in a complex way) on the correlator (2^2) |2C[ |. B is a scaling function discussed in appendix [B|, 
and q is the "plateau value" of B(t), equal to in the zero temperature limit. Meanwhile, the fluctuation dissipation 
theorem, obeyed for t < tf, is violated around t ~ tf , and becomes : 

dB{t)/dt = 2TR{t); t > tf. (2.15) 

The effective temperature T and the plateau va lue q are identical to those obtained in the aging case . 
The velocity- force characteristics is given by (|2.9| ) , and in the limit of small velocities, 



v - F/t b {v). (2.16) 

The time tb(v) plays the role of an effective friction coefficient, controlled by the velocity. The v — F characteristics 
is a power law. 
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(2.17) 



Let us mention for completeness the presence of a third time scale, called t' a in 

consider an exponential correlator, we have exp(— B(t) — v 2 t 2 ) = exp(—B(t)) x exp(— v 2 t 2 ) and in our case t' a ~ v 
We believe that apart from this point, the results of j2(J all qualitatively apply to the exponential correlator, and 
anyway t' a does not play a direct role in the dynamics of the short range correlated models. 



20), defined by B(t' a ) = v 2 t' a . As we 



III. GEOMETRICAL DESCRIPTION OF THE DRIVEN STATIONARY DYNAMICS 



In 1 21 1 was proposed a geometrical description of the relaxational dynamics of the particle. This approach makes 
use of a comoving frame, defined by the eigendirections of the hessian matrix VVV r (x) at the precise point x(£) 
where the particle stands. This frame is made of N vectors {e^}, each one eigenvector of the hessian VVV r (x). The 
distribution of the corresponding eigenvalues is a semi-circle of radius 4, shifted towards the positive values, and such 
that the lowest one is equal to —S. Each eigenvector has an eigenvalue A^ — S, and the density of states of the Xi 
is : 



p(X) = (Stt)- VA (8 -A). (3.1) 

The quantity S is positive, and depends linearly on the energy of the system, e.g. S(t) = 4 + 2V(x(t))/N = 4 + 2£ (t). 
In the aging case, S(t) is a time-dependent function, while in the stationary case, iS is constant. 
One projects the instantaneous velocity x is the above frame such that: 

JV 

x = ^7iei (3.2) 

i=l 

Because the spacing of the eigenvalues is of order 1/N, the set of Xi becomes dense, and one replaces the discrete sum 
over the index i by a continuous one, involving the semi-circular density of eigenvalues p(X). 

,8 

x 2 = / dA p(X) g{X,t) (3.3) 
Jo 

The distribution g(X,t) represents the mean value of the component x 2 , locally averaged over the indices i such that 
Xi ~ A. We have justified in [^jj the following self-similar form for g{X,t): 

x 2 ~ ff(Ai,t) 

g(X,t)=S(t)G(x/S(t)), (3.4) 



where A stands for any direction with a curvature of the potential equal to A — S pi) . The prefactor S in front of the 
distribution comes in fact from an assumption about the value of the exponent — k governing the power law decay of 
£(t) and S(t) ~ t~ K , which we believe to be —2/3. This assumption is supported by our numerical results. 
The characteristic times tf and % are controlled by S, and scale like : 

t b ~S-V 2 ; (3.5) 

making the instantaneous velocity equal to : 

ii(t) = F-S 3/2 = F/t b (t). (3.6) 

When a constant force F is applied, the dynamics changes from an aging linear-response behaviour to a stationary 
regime pl[ . At short times, the displacement u{t) is proportional to the integrated response F ■ TZ{t, 0), while at long 
times, it becomes equal to v ■ t. 

One expects the stationary regime to take over the aging regime when the dynamics is dominated by the external 
force T rather than by the gradient -V^(x(t)). This happens at a time tp, inversely related to the magnitude of 
the force. Our numerical definition of tp, is the time where the slope of the asymptotic curve u(t)/F ~ vt/F is equal 
to the slope of the (logarithmic) integrated response TZ(t, 0), as shown on Figure (111) with F = 0.3. 
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In our situation, the dynamics is controlled by the value of S; inverse friction and diffusion coefficients are both 
proportional to S 3 ^ 2 . In the aging case, S(t) tends to zero as a power law, and the dynamics of the system is slower 
and slower. A look at Figure (||) however shows that in the presence of a force, S does not go to zero, but to a 
finite value S(F), controlled by F, and inversely related to the magnitude of F. The same is true for the energy 
£(t) = — 2 + S(t)/2, which stands higher than in the absence of driving force. Both diffusivity and mobility are kept 
finite thanks to a non-zero driving force. What is needed is to compute S(F). For this purpose, one assumes that the 
self-similar form ([TJ) is still valid in the stationary regime. A justification is provided in appendix |aJ. 

The zero temperature relaxation equation is: 

*i = -diV(x) + T u (3-7) 

while the energy obeys : 

i 

£{t) = — x 2 (t)/7V + F ■ u(t). (3.8) 

One uses now the distribution g(X,t) of the instantaneous velocity components xf, the density of eigenvalues p(X) 
and finds : 

- £(t) + F ■ u{t) = I dX p(X) g(X,t). (3.9) 



In the stationary regime, £{t) — and u(t) = v. The equation (3.E) reduces to a balance between the mechanical 
power given by the force, and a kind of intrinsic dissipation (x 2 ). 



Fv=l dX P (X)g(X). (3.10) 

Assuming that g is still equal to SG(X/S) (cf appendix [A|), one gets : 

F-v~S 5/2 . (3.11) 



From (3.11) and (3.6), one finally finds S as a function of the force, 

S~F 2 , (3.12) 

the resulting velocity force characteristics, 

w-F 4 , (3.13) 

and the force and velocity dependence of the time scales : 

t f ~ F- 2 ~ v~ 1/2 ; (3.14) 
t b ~ F~ 3 ~ v- 3 ' A . (3.15) 

These results are in full qualitative agreement with the findings of Horner [^0) . The main relaxation time tb does not 
scale as as could be expected from a simple dimensional analysis, but is shorter, such that lim„^o v • tb(v) = 0. 

One determines the cross-over time tp by a matching argument. In the linear response regime, S(t) decreases 
as t~ 2 / 3 , as the force acts only as a weak perturbation. The linear response breaks down when the perturbation 
modifies the nature of the relaxation itself. This happens when S(tp) reaches the order of magnitude of its limit value 



S(F) = F 2 (equation |3.12| ), leading to t p 2/3 = F 2 , or 



t F =F- 3 . (3.16) 

Physically, this means that a typical coordinate /j of the force T, along a downhill direction i, is of the same order of 
magnitude than the gradient of the potential —diV , or the instantaneous velocity ij EM. From equation (|3.4[) and 



fi ~ F, one gets F 2 ~ f 2 ~ x? ~ <S, in agreement with (3.12) 



Let us mention that a qualitatively similar cross-over has been observed in the simulated dynamics of a driven 
polymer, in the presence of quenched disorder pgj]. 
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IV. NUMERICAL RESULTS 



We present numerical results which support the findings of the previous section. 

Figure (^) shows S(F) versus F, in log coordinates, and in regular coordinates (inset), for F = 
0.05,0.1,0.2,0.3,0.4,0.5 and F — 0.6. The squares are the values obtained with t max — 200 (h = 0.1), and the 
full curve with t max = 400 (h = 0.2). One sees that the three first values are not well converged. If we except them, 
the overall shape of the curve is concave (downward curvature). The slope of the tangent curve between the arrows 
gives an exponent equal to 1.81 (h = 0.1) and 1.87 (h = 0.2). Because the curve is concave, we believe that these 
values are a lower bound for the real exponent, compatible with our prediction 2. 

Figure (^) shows v(F) versus F, in log coordinates, and in regular coordinates (inset). The squares are the values 
obtained with t max — 200 (h = 0.1), and the crosses with t max — 400 (h = 0.2). As for Figure (||), the three first 
values are not well converged. If we except them, the overall shape of the curve is again concave. The slope of the 
tangent curve between the arrows gives an exponent equal to 3.73 (h = 0.1) and 3.82 (h = 0.2). Repeating the above 
argument, these values are a lower bound for the real exponent, compatible with 4. 

A plot of tp vs F is reported on Figure (|J). Again the system has not reached its asymptotic regime as far as the 
three first values F < 0.2 are concerned. This can be checked by looking at the first derivative u(t) which must be 
constant when t reaches the upper limit of the time window, here t = 400. The fitted value on the straight part of the 
graph, in logarithmic coordinates, gives an exponent —2.72 instead of —3. Again the true asymptotic limit F — > is 
out of reach, due to our limited computer facilities. 

These numerical results are not good enough to prove the exactness of the equations ( 3.12 ) and (|3.13| ). However 
they provide lower bounds which constraint the exponents to be larger than 1.8 for S, and larger th an 3 .8 for v. On 
the other hand, if we assume that we are close enough to the asymptotic regime where equations ( 3.12 ) and ( 3.13| ) 
apply, one expects the real exponents to be not too much different from the above numerical values. In this respect, 
we think that the numerics is in agreement with our findings. As far as the cross-over time is concerned, the numerical 
exponent is —2.72 instead of —3. A larger time window would certainly improve the agreement. 



V. A UNIFIED DESCRIPTION OF THE OUT-OF-EQUILIBRIUM REGIMES 



In the isolated aging regime, at zero temperature, the correlation function obeys, as a particular case of equation (Bl) 
of appendix |b[ as shown in |l9) : 



b(t,t') 



In 



[Ht)\ 
\h(t')J 



where the parametrisation function is related to the time-scale tb by 

t b (t) = h{t)/h'{t) 
As the time scale % is proportional to 5~ 3 / 2 , we have : 

h(t) 



(5.1) 



(5.2) 



exp 



C / ds S 3 ' 2 



h(t>) 

b(t,t') = C I ds S 3/2 



On the other hand, equation (3.6) leads immediately to 



u{t) - u(t') 



= C ds S 3 > 2 . 



(5.3) 



(5.4) 



Now, one observes that the scaling form ( 2.12 ) resemble to ( 5.1 ) , (pl| ), with q — 0. We prove in appendix [Bj tha t 
the scaling function of the aging regime (l{J and the driven regime |2d|l are indeed equal, and thus B(x) = x in (|2.12| ) 
(strictly speaking, B(x) is only proportional to x, but one can choose % such as B(x) = x). The equations ( |5.3|) and 
( |5.4D make sense in the aging regime as well as in the stationary regime. 

The integral J., ds <S 3 / 2 is the effective time variable for the system, interpolating smoothly between ln(£) (aging, 
linear response regime) and t/tb(F) (stationary regime) while S~ 3 ^ 2 is an effective age, growing like the waiting time, 



6 



in the aging regime, and bounded in the stationary regime. Interestingly, a similar effective age has been used in the 
context of the stick-sli p m otion in d ry friction experiments |^| . 

The prediction for ( |5.3| ) and ( |5.4[ ) is checked by plotting b(t,t') vs [u(t) — u(t')]/F, shown on Figures (||) and (Q). 
One expec ts b(t, t' ) and [u(t) — u(t')]/F to be proportional, both in the aging and stationary regimes, provided 
equations (5Ji_5.4) hold, which is the case for a time separation t — t' large enough. 

On Figure (|6), b(t,t') is plotted against [u(t) — u(t')]/F for F = 0.1 (crosses, squares and diamonds) in the linear 
response regime and for F = 0.5 (continuous lines) in the non-linear regime. The force is zero till t = t' , and then 
switched on; t' takes the value 0, 20 and 40. As far as F — 0.5 is concerned, the transition from linear to non- 
linear regime is not visible on this curve, and in any c ase very smo oth. The slope of the curve defines the effective 
temperature 2T, equal to the ratio C/C in equations (5.3) and (5.4). The effective temperature thus makes sense in 
both linear and non-linear regimes. 

As the force is switched on at t', there is a short-time "elastic" displacement. This is how the directions with a 
positive curvature respond to the new static constraint, and this corresponds to the short horizontal step at the origin, 
seen on Figures (inset) and (0). The finite slope part of the curv e cor respo nds to the slow wandering motion of 
the particle in the energy landscape, in the regime where equations (5J3) and (5^) apply. Thus, we conclude that 
Figure (^) support the proportionality of u(t) — u(t') and b(t,t'), once the short time regime has been taken into 
account. 

A close look near the origin of the graph (inset of Figure (^) shows that the F — 0.5 curve is slightly shifted from 
F = 0.1, but parallel to it. This shift goes rapidly to zero as F — ► 0. The shift is presumably there because 0.5 is 
already a large value of the force, leading to a departure from the ideal curve corresponding to F <C 1. 

Figure (0) is the same as Figure @ for F = 0.01, F = 0.1 and F — 0.5, for three values of t', 0, 20 and 40, and gives 
additional details on the short time response of the particle. Again, the horizontal part of the curves corresponds to 
the short-time displacement ("elastic" or reversible) while the finite slope regime corresponds to the slow motion in 
the energy landscape ( "plastic" or irreversible) . 



VI. CONCLUSION 



In this paper, we have proposed a consistent picture for the stationary driven dynamics, in the mean field approxi- 
mation and zero temperature limit, of a particle in a quenched, exponentially correlated, random potential. 

The velocity (v)- force (F) relation is a power law v — F 4 , while the main relaxation time scales as tb — v~ 3 ^ 4 . 
The product v ■ t\, tends to zero as v vanishes. These findings are consistent with earlier work |po| . The driving force 
is found to generate a relaxation time smaller than the "dimensional" time scale v^ 1 , which is probably a generic 
feature of the mean-field short-range correlated potentials. 

If the force F is small enough, a linear response around the aging regime is found, up to a time ip, scaling as F~ 3 . 
A plot of the displacement (u(t) — u(t'))/F vs the correlation b(t, t') shows no sign of discontinuity, when the linear 
response regime is replaced with the non-linear stationary regime. We interpret it by saying that, when a small force 
is applied, the dynamical properties of the system (mobility, diffusivity) are controlled by the effective age S~ 3 ^ 2 . The 
quantity 6> 3 / 2 is proportional to the number of negative eigenvalues in the spectrum of the hessian of the hamiltonian. 
The effective age is proportional to the waiting time in the aging regime, and finite in the stationary case. 

The effective time J* S 3 / 2 ds, closely related to the correlation function b(t,to), grows logarithmically with t in the 
aging regime, and linearly with t in the stationary regime. The effective temperature T generalising the fluctuation 
dissipation theorem, remains unchanged in the non-linear regime. However, the geometrical meaning of T, if any, is 
still unknown. 

Future work will determine to what extent are the present features generic from other short range correlated models, 
and finite dimensional models. Even though such a power law dependence of the characteristic times in the driving 
force is not observed in realistic systems, the qualitative behaviour presented in this study -cross-over between linear to 
non-linear regime, coexistence of aging and non-linear stationary dynamics-, could indeed be a very generic situation. 
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APPENDIX A: THE ENERGY BALANCE 



Let Yi(i) be the coordinates of the instantaneous velocity x(i) in the comoving frame {ej(i)} (pT^). When T = 0, 
this definition is equivalent to say that 7^ is the coordinate of — VU. Using the local average defined in [|IJ, one finds: 



|x|| 2 (tH£x^)=£ 7l 2 W, 

i i 

= N fd\ p(X) g(\,t). 



The derivative of ||x|| 2 (i) reads : 



N 



d t \\±f(t)=Y,d t (-d 3 V + T 3 f 

3 

= -2^2d jk V -xj -x fe 

jk 



= -2N J dX p(X) (X-S)-g(X,t). 
We deduce that, in the stationary situation, for all S, 

dX p(X) X .9(A) = S [ dX p(X) <?(A), 



which is in favour of a scaling form g(X) = TG(X/S). 
As diV(x(t)) = —Xi(t) + Ti, the equation for £(t) is : 



£(t) = N- 1 J2d J V. Xj 

3 

= N- 1 J2{-^j(t)+^-^3<t)} 



The product N 1 J^j J~j ' *j 1S by construction equal to F ■ u(t). Thus, (this is equation 3.9): 



£(t) 



dX p(X) g(X, t) + F ■ u(t). 



The energy balance ( 3.11 ), and the factorised form of g(X) imply in the stationary regime : 

T • S 3/2 oc F 2 ■ S 3/2 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



However the relation between S and T remains undetermined by the present argument. For the sake of simplicity, 
we can suppose that the equality S = T, true in the aging regime, remains true in the stationary regime. This 
assumption is in fact equivalent to a matching argument, when the distribution g(X,t) = S(t)G(X/S(t)), crosses over 
the distribution g(X) = F 2 G(X/S(F)) around t = tp. The matching of g(X,t) and g(X) leads to the identification 
S(F) = F 2 . One cannot rule out, rigorously, more complicated behaviours, which could lead to a different velocity- 
force characteristics. The assumption r = S is just the most natural one. 



APPENDIX B: THE SCALING FORM OF THE CORRELATION FUNCTION 



In the isolated situation, the correlation function in the aging regime reads, for any finite temperature T |19 3 



b(t,t') = q + B 



. \Ht')J 



(Bl) 
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An general equation for B(u) is obtained in reference |i9[ | (equation 6.22, with the opposite sign convention for /), 
and reads: 

= B(u)f"(q)-f'(q + B(u))+f'(q) 

+ ^f"(q) £ du' B'(u')f"(q + B{u))B{u - u>) (B2) 

whose solution is B{u) — C st xu, leading to (5T), with q = 0. The function / = exp(— x) stands for the correlator (|2.2|). 
T for the temperature, q for plateau value of the correlation function 6, and x f° r the fluctuation dissipation violation 
parameter (see equation B4 below). 

On the other hand, in the stationary regime, the equation for b(t, t') = B(t — t') = q + B(t — t') and R(t—t') = r(t, t') 
is (equation (2.9) in reference (^0|, again with the opposite sign convention for /): 

d t B{t) = 2T - Qf ds Af"{B(s) + v 2 s 2 )R(s^j ■ B(t) + J ds Af"(B(s) + v 2 s 2 )R(s)B{t - s) 
+ ds { (Af'(B(t + s) + v 2 (t + s) 2 ) - 4f'(B(s) + v 2 s 2 )^R(s) 

+ (if"(B(t + s) + v 2 (t + s) 2 )R{t + s)- if"(B(s) + v 2 S 2 )R{ S )^)B{ S )} (B3) 

One knows that the main relaxation scale tb is much smaller than v^ 1 , and asymptotically, lim t _^o v ■ % = 0. The 
above integrals can be safely cut beyond a cut-off A such that % <C A <c v^ 1 . The contributions J. are negligible 

because the relaxation of B(t) has already taken place, while in the integrals J^, the term v 2 s 2 can be neglected 
compared with B(s) in the argument of the correlators /' and /". 

One introduces the quasi fluctuation dissipation parameter X , defined by: 

- X(B(t)) ■ dB(t) /dt = R(t). (B4) 

X is equal to its equilibrium value — 1/2T if B < and to \ if B > 0. Equation ([BO]) becomes: 

d t B(t)=2T+ ^ ds4:f\q + B(s))X(B(s))dB(s)/dsJB(t)-J^ 4f"(B(s))X(B(s))d~B(s)/ds ■ B(t - s) 
- £ ds { (4f'(q + B(t + s)) - if'(q + B(s))) X(B(s))dB(s)/ds 

+ Uf"{q + B(t + s))X(B(t + s))dB(t + a) /6s - Af"{q + B(s))X(B(s))dB(s)/ds) }s(t) (B5) 

Each integral J has to be split to take into account the short time quasi-equilibrium regime and the long time regime. 
As the time scale tf separates these two regimes, one writes J Q = J a a+t/ + J a+ * f + Jb-t f - ^ ne P aram eter X is then 
set to — 1/2T or \ accordingly, and most of the integrals can be reduced to boundary terms. One neglects the time 
derivative dtB(t) in the asymptotic long-time regime, and the result is: 



= 2T+(^ X J t q f"(q + B)dBj(q + B(t)) 
~ 4q ( X+ 2^) X + £(*))) 



-A X f ds dB(s)/ds f"(q + B(s))(q + B{t-s)) (B6) 
By using lim t _ 00 !B(t) = oo, q 2 f"{q) = T 2 and -AxJf"(B)dB = 2T/q, the equation (^|) for B coincides exactly 



with (B2). As the equation (B6) is invariant upon time dilatations, B(u) = B(t/tb) is a solution of (B2) and without 
loss of generality, one has: 

B(u) = B{u) = u, (B7) 

which is the announced result. 
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Captions 



FIGURE 1. Determination of the cross-over time tp, denned as the time where the slope of the integrated response 
lZ(t, 0) is equal to the velocity 1/Flim^oo ii(t). 

FIGURE 2. Family of curves £ (t) + 2 for increasing forces, ranging from F = 0.05 to 0.5. The limit value 
lim^oo £ (t) + 2 is a monotonically increasing function of F, equal to S(F)/2. The effective mobility and diffusivity 
are directly related to S(F). The system stays above the marginal states, in a region with a finite extensive number 
of downhill directions. 

FIGURE 3 The parameter S as a function of the force, for F — 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6, in log coordinates, and 
normal coordinates (inset). The boxes stand for a run, up to a time t — 200 while the straight line corresponds to 
t = 400. Whenever the boxes differ from the line, the value is not converged. See text for details. 

FIGURE 4. The velocity v as a function of the force F, in log coordinates, and normal coordinates (inset). Same 
remark as for Figure (4). 

FIGURE 5. The time tp as a function of the force F in log coordinates. Inset : tp as a function of F. The three 
first values are not accurate (tp larger than our maximum time). The fitted exponent of the straight part is —2.72 
instead of 3; —2.72 is a lower bound for the real value. 

FIGURE 6. The correlation b(t, t') vs the displacement [u(t) - u(t')]/F, for F = 0.1 and 0.5. The force is switched 
on at t', successively equal to 0,20 and 40. Inset: the short-time behaviour. See text for details. 

FIGURE 7. Same as Figure (6), with F = 0.01, F = 0.1 and F = 0.5. 
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FIG. 1. 
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[u(t)-u(f)]/F 

FIG. 6. 
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